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Tunneling of a quantum particle to the classically inaccessible region 

J---. ' in the intrinsically degenerated system is investigated here by means of 

quasienergy eigenstates. The exact resonance Su> = u> — co c = and near res- 

O ■ onance 5uj ^ cases are explored numerically. It is shown that in both cases 
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all quantum states are localized. Correspondence of the quantum dynamical 
barriers to the classical invariant curves is demonstrated. The phenomena 
considered in this article can be observed when an ultrasound wave propa- 



^^ ' gates perpendicularly to a magnetic field and interacts with a 2D electron gas 

(N 

r™ in a semiconductor heterostructure. 

In this paper we treat the role of quantum dynamical barriers in localization of quantum 
P-i. 

states at the cyclotron resonance. Dynamical barriers in quantum systems are counterparts 

& . 

of classical invariant tori (or KAM tori) which block diffusion of a particle along the phase 

space. Similarly, in quantum case dynamical barriers inhibit the diffusion of a wave packet 
in a Hilbert space that leads to localization of quantum states. This phenomenon was 
previously analysed for the kicked rotator model [1, 2] and the nonlinear oscillator in a 
monochromatic radiation field [3]. However quantum dynamics was explored maily for the 
systems, whose classical counterparts are accidental degenerate and the conditions of the 
KAM theorem are satisfied. The dynamics of the intrinsic degenerate system (kicked simple 
harmonic oscillator) was explored by Berman, Rubaev, Zaslavsky [5] under the condition of 
strong and weak chaos of the classical limit. In this work by using quasienergy eigenstates 
we treat tunneling through the dynamical barriers and localization of quantum states in an 
intrinsically degenerated system where the KAM theorem is invalid. 

Our system is a charged particle subject to both a uniform magnetic field and a field 
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of a longitudinal monochromatic wave, propagating perpendicularly to the magnetic field 
direction. The phase space of the classical system at the condition of the cyclotron resonance 
uj = uj c , where u and uo c are the wave and cyclotron frequencies, consists of the infinit 
number of cells comprising closed classical orbits and separated by a separatrix lattice [4]. 
In the classical case in the resonance approximation a particle has no possibility to penetrate 
through a separatrix from one cell to another. In this paper we will show that in the quantum 
case the dynamical barriers become transparent for the wave packet even in the resonance 
approximation. 

The theory of nonlinear resonance in the intrinsically degenerated system was developed 
by G.A.Luna-Acosta and the authors in Ref. [6]. It was shown that the Hilbert space of 
the discussed system breaks up into quantum resonance cells. The boundaries of these cells 
in the quasiclassical limit correspond to the the separatrices in the classical phase space. In 
this paper we numerically investigate the role of these "quantum separatrices" in localization 
and tunneling phenomena at the cyclotron resonance. 

The Hamiltonian of the considered system has the form 

H = VP c ; +V Q cos(kx - out) =H + V(x, t), (1) 

where e and m are the charge and mass of a particle, c is the light velocity, p is the 
momentum, k and Vq are the wave vector and amplitude of the perturbation. It is convenient 
to choose the gauge of A in the form A = (0,Hx,0) in order to produce the magnetic field 
H along the z-direction. In this gauge the momentum p y is an integral of motion, hence we 
have to determine the dependence of the wave function only on two variables, x and t. 

Since the perturbation is periodic in time, the Floquet theory can be used for description 
of time evolution of the system in terms of quasienergy (QE) spectra and QE eigenfunctions 
ip q {x,t). The QE eigenfunctions are eigenstates of the evolution operator for one period of 
oscillation of the external field T = 2-k/uo. They can be defined by (see for example Ref. [7]) 

iE t 
il> q (x, t) = exp( -2-) J2 Yl Aq n,s^n{x) exp(-isut), (2) 
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where E q is the QE eigenvalues, ip n {x) is the nth eigenfunction of a simple harmonic os- 
cillator, h is the Plank's constant. After substituting Eq. (2) into the time-dependent 
Schrodinger equation, we obtain the system of the uniform algebraic equations: 

(E q - hu c n + kuJs)A q ns = J2 Vn,m( A m,s+l + A m,s-l), ( 3 ) 

m 

where the quasienergy E q is counted from the ground state huj c /2, the matrix elements V n ^ m 
was defined in Ref. [6]. 

First let us consider the exact resonance case where 5uu = uu c — uj = 0. If Vq = each 
definite value n — s in Eq. (3) corresponds to an infinitely degenerated QE level and these 
levels are separated by the distance hu. When Vq ^ the degeneracy is broken and each 
initially degenerated level splits into QE spectrum. All the spectra are identical, and the 
distance between the extreme levels in one spectrum has the order of Vq. If Vq <C Tiuo we 
can, at first order, neglect the interaction between the individual spectra and treat one of 
them. 

QE eigenvalues E q and eigenfunctions A q n s were calculated in Ref. [6] for resonant 
(uj = uj c ) and near resonant (5uj ^ 0) case by the perturbation theory assuming VqJTiuj << 1 
and h5ujn <^C hu for all studied values of n. In the first order approximation Eq. (3) yields 

(E q - h5uon)A q n = V n , n+1 A q n+1 + K,n-i^-i, (4) 

where A q n = A q in . The matrix elements V Ujn ±i are symmetric, real and for n ^> 1 can be 
expressed via the Bessel functions of the first order 
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K,n+i = yW^— e iJiiV^nh), (5) 

where h = (ka) 2 , a = (|§) 1//2 is the magnetic length. When 5u = 0, the amplitude of 
the wave Vq, which is measured in units of hu, does not influence the eigenvectors A q n but 
changes only the energy scale, hence the effective Plank's constant h is the only parameter 
in this problem. The oscillating matrix elements V nyn+ i define the quantum resonant cells 
- relatively independent dynamical regions in the Hilbert space [6]. Each quantum cell 
corresponds to two classical cells (in action-angle variables). 



When 5uj 7^ and Vq = according to Eq. (4) the levels are separated from each other 
by the distance h5ui. If Vq 7^ interaction between the levels appears, however, due to Eq. 
(4) this interaction is essential only within the region nhbio ~ Vq, whereas the levels satisfied 
the condition nhbuj ^> Vq practically are not affected by the perturbation. 

Let us discuss the main features of the QE spectrum which substantially define dynamics 

of the system. In the exact resonance case the spectrum E q is almost equidistant near the 

top and the bottom at the distance hCo [6]. The frequency Co in the quasiclassical limit is 

identical to the frequency of small oscillations near the center of the resonance in the phase 

space. In order to derive Co for the near resonance case, let us expand f(n) = V njn+ i near 

some point n , defined below, into Taylor series and introduce the differential operator by 

d 2 
A n+ i + A n _ 1 = -—A n + 2A n . (6) 

Then the finite-difference Shrodinger equation (4) takes the form 

[£ q - (n - n )(HSoo + 2f(n )) - 2/(n )] A$ = /(«o)^^ + f"(n )(n - n ) 2 ^, (7) 

where £ g = E q — h5um . If the condition 

hSco = -2/'(no) (8) 

is satisfied equation (7) takes the simple harmonic oscillator form with the center of oscilla- 
tions no defined by Eq. (8). According to Eq. (8), when 5oo 7^ there is only a finite number 
of the resonances. If Vq is small, i.e. Vq < h5uj, the unperturbed levels nhSuo are not coupled 
by the perturbation (see Eqs. (4) or (8)) and the resonance disappears. The same is true 
for the classical case. In the limit h — > 0, n — ► 00, and finite nh = I, the positions of the 
resonances and the frequency, defined by Eqs. (8) and (7), correspond to the fixed points 
and frequency of small oscillations near the center of the resonance in the phase space [8]. 

An arbitrary solution of the time-dependent Schrodinger equation with the Hamiltonian 
(1) can be expanded in the simple harmonic oscillator basis 

i/)(x,t) = ^2C n (t)if; n (x)exp(-iE n t/h), (9) 



were E n is the energy of nth Landau level. Coefficients C n (t) determine the time devel- 
opment of the probability distribution defined by |C n (t)| 2 . Instead of integration of the 
time-dependent Schrodinger equation we express the coefficients C n (t) in terms of the Green 
function G njTl i 

C n {t) = } J G n y (t, t )C n > (t ) , (10) 

n' 

which can be determined by using the complete set of the QE eigenstates 

G n , n .(t,t ) =Y. Aq n A l>Z- iEq{t - t0)/h - (11) 

Thus, we have the following procedure for computing the probability distribution at time 
t which is used below. First, by numerical solving the finite-difference Schrodinger equation 
(4) we compute the QE eigenfunctions A q n and the QE eigenvalues E q , then using the QE 
eigenstates we compute the Green function (11), and finally we use the Green function for 
calculating coefficients C n (t). Since the dependence of the coefficients C n (t) on time t in 
Eqs. (10), (11) is explicit, they can be easily computed at arbitrary time t. 

We have performed the numerical experiments in order to explore how the dynamical 
barriers in the intrinsically degenerated system (web tori [9]) influence the evolution of 
quantum states with various initial conditions at the cyclotron resonance. The system 
consists of the weakly interacting quantum resonant cells [6] defined by the matrix elements 
V n ,n+i- Our purpose is to treat electron tunneling between the cells. 

As mentioned above, the QE states completely define the dynamics of the system. If the 
initial state is any QE eigenstate, i.e. C n '(t ) = A q n , , then in accordance with Eqs. (10) and 
(11) the probability distribution does not depend on time t. In particular, if the initial wave 
function A q n , is localized in one cell then it will be localized in the initial cell for any time 
t. The localized QE eigenfunctions correspond to the eigenvalues near the top and bottom 
of the QE spectrum. By employing a coherent state representation (Husimi function) we 
can show that the localization of a wave packet in the Hilbert space entails the localization 
of the Husimi function in the phase space. The contour plots of two Husimi functions 



corresponding to the extream eigenvalues from the top and bottom of the QE spectrum are 
plotted in Fig. 1 (a) (the upper part of the figure corresponds to the top eigenvalue and the 
lower part corresponds to the bottom one). The symmetry of these Husimi functions can 
be explained from Eq. (4) which is invariant under the transformation 

E q -> -E q , Al -> (-l) n Al (12) 

corresponding to the transformation x — ► —x (or out — ► out + n) in Eq. (2). The contour 
plots of the Husimi functions correspond to the classical orbits in two cells of the phase 
space which are shown in Fig. 1 (b). Another kind of initial state is an excited Landau 
level C„(0) = S ntno . In this situation, the initial wave packet contains all the QE eigenstates 
and one can expect that such a state will be extending most intensively, tunneling to other 
resonant cells. 

In our calculations the system involved 432 Landau levels which formed 7 resonance 
cells in the Hilbert space, shown in the upper part of Fig. 2. The initial state in the form 
Cn(0) = 5 n ,n w as placed in the center of the first cell. After the time T ~ 2-k/Co the wave 
packet spreaded over the initial cell. In parallel with the fast dynamics within the initial cell 
a slow process of propagating the probability distribution to the subsequent cells took place. 
The characteristic time of the probability distribution spreading over all the considered cells 
has the order of 27i/u m i n , where w m i n is the minimal distance between the QE eigenvalues 
in the QE spectrum calculated for seven cells. Two snapshots of evolution of the probability 
distribution \C n \ 2 (in a logarithmic scale) as a function of the Landau number n are shown in 
Figs. 2 (a), (b). It is seen that the probability distribution in average decreases exponentially, 
but sharp decay occures only at the boudaries of the cells that verifies our assumption of 
the role of these boudaries as dynamical barriers to the probability flow. 

The probability distribution successively penetrates through the dynamical barriers from 
one cell to the subsequent ones, travelling along the Hilbert space. After filling some cell 
the wave packet reflects from the barrier and interferes (in the sixth cell in Fig. 2 at the 
moment (a) t — 4 * 10 5 and (b) t — 7 * 10 5 , time is measured in units of T = 2tt/uj). 
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As a result, oscillations of the probability distribution appear. These oscillations can be 
explained also from another point of view. Let us consider tunneling of some initial wave 
packet, concentrated on a level with Landau number n in an initial cell, to another level 
with number n', situated in another cell. The main contribution to the Green function (11) 
is given by the QE eigenf unctions with those numbers q, for which both A q n and A q n , are 
sufficiently large or, in other words, any of these QE eigenfunctions must occupy two (or 
more) resonant cells. Such eigenvectors correspond to the eigenvalues near the center of 
the QE spectrum, i.e. the point E q = [6]. As follows from our calculations, the most 
delocalized QE eigenfunctions at the same time turn out to be most oscillating (by module) 
functions of n. The probability distribution in Fig. 2 also oscillates with the minimal period. 
It is interesting that the probability distribution in average is greater near the boundaries of 
the cells (by 1.5-3 orders in Fig. 2). Since these boundaries correspond to the separatrices 
[6] (web tori), this quantum phenomenon is akin to the diffusion of a classical particle along 
the separatrice lattice within exponentially small stochastic regions in the phase space when 
nonresonant terms in the classical Hamiltonian are taken into consideration [9]. 

The probability distribution in Fig. 2 is approximately the same within a cell, hence 
it is reasonable to consider average probabilities, defined by -Pj(t) = J2n t \C ni (t)\ 2 , where n» 
takes the values within the ith resonant cell. The evolution of P\{t) and P2CO (i n the first 
and second cells) is presented in Fig. 3 (the time axis is plotted in a logarithmic scale). 
For the calculations we used the system with the same parameters as in Fig. 2 involving 
seven resonance cells. We suppose that the first two cells do not " feel" the limited boundary 
conditions because the probability distribution at the boundary of the system (the seventh 
cell) is exponentially small and including of more number of cells does not affect Pi (t) and 
P 2 (^)- The probability distribution in Fig. 3 does not leak out to the subsequent cells 
in the limit t — ► 00, hence we can conclude that the quantum states are localized. The 
localization of the quantum states suggests a conjecture (Ref. [6]) about discreteness of 
the QE spectrum for the unbounded system when n — »■ 00. Time- averaged probabilities 



Pi(t) versus cell number for two values of the effective Plank's constant h are presented 



in Fig. 4. In the logarithmic scale the exponential behaviour of these functions is fairly 



evident. As one can see in Fig. 4 the relation P 2 (t)/Pi(t) has the order of 1CT 2 . It should 
be stressed that the sharp decay of the probability distribution occures not everywhere, but 
only at the boundaries of the quantum resonance cells. Thus, we deal with the new kind of 
localization, namely, localization over the quantum resonant cells. The localization length in 
the resonance approximation does not depend on the amplitude of the wave Vq. It is defined 
by the only parameter h = (ka) 2 which determines the behaviour of the Bessel function 
in the matrix elements in Eq. (5). Our numerical experiments have shown that slope of 
the curve on Fig. 4, defining the localization length, in general increases with increasing h. 
However, the decay of the probability distribution at different boundaries between the cells 
changes randomly. The nature of this randomness is discussed below. 

Next effect, which we consider here, follows from discreteness of the quantum number n 
which labels Landau states in the Hilbert space. If we choose the parameter h so that 



V no , no+1 ~ My/^oh) = (13) 

then in accordance with Eq. (4) transitions through the level Uq are blocked and the Hilbert 
space is devided into two disjoint parts. If we will change the parameter h, the condition 
(13) will be successively satisfied for the n + 1-th, n + 2-th, . . . levels. Thus the tunneling 
probability through the dynamical barrier goes periodically to zero. Let us define the pen- 
etration coefficient by P = J2 n >n \Cn{T)\ 2 , where n is located at the boundary between 
the second and the third cells, and place the initial state in the form C n (0) = 5 n y at the 
center of the second cell . The oscillations of P with changing the parameter 1/h ~ H 
are shown in Fig. 5, where the above-mentioned zeroth values of the penetration coefficient 
have been replaced by small numbers. The oscillations have the period AH = 2hck 2 /eb 2 , 
where b is a root (the second one in our case) of the Bessel function (5). The amplitude of 
the oscillations decreases with decreasing h, that does not contradict to the quasiclassical 
formula P ~ exp(— oc/h). The definition of more explicit form of the dependence of the 
tunneling probability P on h is complicated by the wild oscillations of the dynamical barrier 
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permeabilities as a function of h. 

The dynamical barriers become impermeable for definit values of h only in the reso- 
nance approximation when the condition Co <C uo is satisfied and we take into account only 
transitions between the nearest levels. 

Next we present the results of cumulations for the near resonance case Su 7^ which 
are illustrated in Fig. 6 where the probability distribution spreading after the time t = 10 6 
is shown. There is a significant difference between the near resnance and exact resonance 
cases. Though the matrix elements (5) are the same as in the exact resonance case, the 
additional term nhSui in (4) deforms the resonance cells with small n's and destroys the cells 
with large n's. That is indicated in Fig. 6 where the boundary of the probability distribution 
spreading, pointed out by the sharp decrease of this quantity, does not correspond to the 
boundary of the second cell in the case 5u = 0, marked by the second arrow on the figure. 
As follows from calculations with decreasing Vq (or increasing 5u) the boundary of the 
probability distribution spreading decreases. The same is characteristic for the classical 
acidentally degenerate system [8] where in the case 5uo 7^ the resonance occures only in 
several resonance cells in the phase space and the resonances are separated from each other 
by invariant curves. According to equation (8), there is only one resonance at choosen 
parameters, and the probability distribution does not spread via the cells, no matter how 
long we observe the dynamics. This is natural because, as was shown above, the wave with 
a small amplitude Vo does not affect the states with large Landau numbers n. That can be 
also understood from the expression for the Green function (11). As one can see from Eq. 
(4), all the QE eigenstates with Landau numbers n from the region nh5u ^> Vq are localized, 
A\ = 5 n)q , with the characteristic localization length An ~ Vo/h5uj (that can be shown for 
example in the quasiclassical limit). On the other hand in order for any state to evolve from 
level n to level n' both A q n and A q n , must be large, but if A q n = 5 nyq , then A q n , = for all 
n' 7^ n. Hence boundary of the probability distribution spreading should be determined from 
the condition of localization of the QE eigenfunctions (the condition nhbuo > Vq was satisfied 
for the boundary of the probability distribution spreading in Fig. 6). The steep decay of 



the probability distribution in Fig. 6 in the quasiclassical approximation corresponds to 
absence of tunneling between invariant curves. A qualitative comparison of permeabilities 
of the invariant curves in the case 5u 7^ with permeabilities of the "quantum separatrices" 
in the case 5uj = (see Fig. 2) indicates that the latter are much more transparent (compare 
also with Ref. [1]). 

In conclusion we would like to do the following remark. The nonperturbed Hamiltonian 
Ho yields an unlimited (from above) set of Landau levels with equal distance hu c between 
them. When we add a monochromatic wave with the same frequency uj = uj c , a particle 
can pass to other levels and an interesting question arises: is the diffusion limited or not? 
The answer depends on how the perturbation depends on coordinate x. If the dependence 
is linear then quantum states are delocalized and the particle goes to infinity in the limit 
t — > oo (as in the classical case) [10, 11]. As follows from the results of our work, thanks to 
nonlinearity of the potential Vq cos(kx — uot) the quantum states (as classical ones) turn out 
to be localized. However, the peculiarity of localization in the latter case is that the sharp 
decay of the probability distribution occures not everywhere, but only at the boundaries of 
the quantum resonance cells. 

In the next paper we will focus on studying destruction of the quantum resonances 
under the influence of nonresonant terms in the Hamiltonian (1) when the amplitude of 
perturbation is not small. 
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FIG. 1. Contour plots in the coherent state representation for the eigenstates with eigenvalues 
at the top and bottom of the QE spectrum N = 100, h = 0.52 (a). Poincare surfaces of section for 
one cell for the same parameters of the Hamiltonian (b). 

FIG. 2. Probability distribution for seven quantum cells, h = 0.6, Vq = 0.1, the number of levels 
N = 432, (a) t = 4 * 10 5 , (b) t = 7 * 10 5 . The initial state was placed in the center of the first cell. 

FIG. 3. The average probability distribution for (a) the first cell P\(t) and (b) the second cell 
P^t). The parameters and the initial conditions are the same as in Fig. 2. 

FIG. 4. Time-averaged probability distribution for two values of the effective Plank's constant 
h versus cell number. The parameters are the same as in Figs. 2, 3. 

FIG. 5. Oscillations of the penetration coefficient from second to the subsequent cells, t = 4*10 4 , 
iV = 100 (three cells), V = 0.1. 
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FIG. 6. Probability distribution for the near resonance case 8u = 0.003, t = 4 * 10 6 , h = 0.52, 
Vq = 0.1, the number of levels N = 100. The initial state in the form C n (0) = 5 n ^ no was situated 
at the level no = 6. 
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